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Abstract. We analyze the temporal behavior of Network 
Bright Points (NBPs) using a set of data acquired dur- 
ing coordinated observations between ground-based obser- 
vatories (mainly at the NSO/Sacramento Peak) and the 
Michelson Doppler Interferometer onboard SOHO. 

We find that, at any time during the observational 
sequence, all the NBPs visible in the NaD 2 images are 
co-spatial within 1" with locations of enhanced magnetic 
field. The "excess" of NaD 2 intensity in NBPs, i.e. the 
emission over the average value of quiet regions, is directly 
related to the magnetic flux density. This property implies 
that, in analogy with the Ca II K line, the NaD 2 line center 
emission can be used as a proxy for magnetic structures. 

We also compare the oscillation properties of NBPs 
and internetwork areas. At photospheric levels no differ- 
ences between the two structures are found in power spec- 
tra, but analysis of phase and coherence spectra suggests 
the presence of downward propagating waves in the in- 
ternetwork. At chromosphcric levels some differences are 
evident in the power spectrum between NBPs and inter- 
network. At levels contributing to the NaD 2 emission the 
NBPs show a strongly reduced amplitude of oscillations 
at the p - mode frequencies. At levels contributing to the 
Ha core emission, the amplitude of network oscillations is 
higher than the internetwork ones. The power spectrum of 
NBPs at this wavelength shows an important peak at 2.2 
mHz (7 minutes), not present in the internetwork areas. 
Its coherence spectrum with Ha wings shows very low co- 
herence at this frequency, implying that the oscillations at 
these chromospheric levels are not directly coupled with 
those present in lower layers. 

Key words: Sun: photosphere - Sun: chromosphere - 
Sun: magnetic fields - Sun: oscillations 



1. Introduction 

The chromospheric bright network has long been observed 
in narrowband spectroheliograms taken in the H and K 
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cores of the Ca II resonance lines. The Ca II network typ- 
ically shows H and K profiles with high double peaks and 
enhanced line wings that persist for extended periods of 
time (longer than 10 minutes, see, e.g., Rutten & Uiten- 
broek 1991 ). The chromospheric network emission pattern 
is cospatial with small-scale magnetic field concentrations, 
and defines the supergranular network boundaries. It is 
this atmospheric component that produces the correla- 
tion between H and K excess line-core flux and magnetic 
activity of cool stars ( ^chrijver et al. 198E ) . 

The dynamics of the network elements, compared with 
the internetwork or quiet chromosphere, has been exten- 
sively studied (especially from the observational point of 
view) since these small-scale structures can be important 
in channeling the energy from photospheric layers to the 
transition region and corona ([Kneer fc von Uexkiill 1986, 



1993, pcubner fc Fleck 1990| , [Kulaczcwski 1992] , Al et al. 
1998). 

An assessment of the spectral characteristic properties 
of Network Bright Points (NBPs) at different layers in 
the atmosphere has been provided by Lites et al. (1993) 
using spectral observations in the range of the Ca II H 
line. In their work, these authors analyzed spectrographs 
observations of a single network bright patch and of sev- 
eral internetwork points. The wavelength shifts of photo- 
spheric and chromospheric lines allowed them to perform 
a compared analysis between the dynamics of the two at- 
mospheric components. One of the relevant characteristics 
they describe is that at chromospheric levels (Ca II H3) 
the NBPs show long period oscillations (y < 3 mHz) not 
correlated with oscillations in the lower atmosphere, while 
they lack power at higher temporal frequencies. The inter- 
network regions display instead enhanced power at higher 
frequencies, well correlated with photospheric oscillations. 
The presence of these low frequency oscillations in the 
network has been confirmed by Lites (1994) also for the 
chromospheric He I 10830 line, in contrast to Bocchialini 
et al. (1994) which observe, for the same line, oscillations 
only in the 5 minutes range. 

An enhanced power in the low frequency range for net- 
work points with respect to the internetwork has also been 
observed by Kneer & von Uexkiill (1986) in the center of 
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the chromosphcric Ha line. These authors however inter- 
pret this feature as not due to oscillations, but of mainly 
stochastic origin, and attribute it to erratic motions of the 
corresponding photosphcric footpoints. 

The problem is still open, and further observations to 



better address this issue are required (Lites 1994). In par- 
ticular one would need observations: 1- on a larger num- 
ber of NBPs, to improve on the statistics; 2- at different 
heights in the atmosphere, since the analysis of the co- 
herence between fluctuations at different levels can help 
exploring the nature of oscillations. To this end, a reliable 
method for the identification of the same physical struc- 
ture at different atmospheric levels is mandatory, since 
the inclination of the magnetic field could displace the 
chromospheric network points with respect to the corre- 
sponding photospheric ones. 

In this paper, we address some of these issues, and 
present observational results on the NBPs and internet- 
work characteristics as derived from a multiwavelength 
analysis. The observations were obtained in August 1996, 
during a coordinated observing program between ground- 
based observatories and the Solar and Heliospheric Ob- 
servatory (SOHO) . For the ground-based observations we 
used the cluster of instruments at the NSO/Sacramen- 
to Peak R.B. Dunn Solar Telescope (NSO/SP-DST), that 
could provide a complete coverage at lower atmospheric 
levels. The dataset used is described in Sect. |[ General 
properties of a sample of NBPs, followed from the photo- 
sphere up to the chromosphere and including their rela- 
tionship with the magnetic structures, are given in Sect. 
||. The temporal development of the NBPs is described 
in Sect. [|. Sect. || and || provide an analysis of the power, 
phase difference and coherence spectra for the fluctuations 
observed separately within the NBPs and the surrounding 
internetwork. Finally, discussion and conclusions are given 
in Sect. 0. 

2. Observations and data reduction 

A description of the general data acquisition has been 
presented in Cauzzi et al. (1997, 1999). Table 1 gives a 
summary of the observing setup, and we recall here only 
some short information on the data used in this paper. 
Monochromatic intensity images were obtained with the 
tunable Universal Birefringent Filter (UBF) and the Zeiss 
filter, at high spatial and temporal resolution. Several 
spectra were obtained around the chromospheric Call K 
line with the Horizontal Spectrograph (HSG). The spec- 
tra have been acquired setting the spectrograph slit on 
different bright points at different times; the field of view 
in the HSG row of Table 1 hence refers to a single slit 
exposure. Onboard SOHO, the Michelson Doppler Imager 
(MDI, (Bchcrrcr ct al. 1995 ) acquired data in high resolu- 
tion mode, i.e. with an image scale of 0.605"/pixel. Maps 
of pseudo-continuum intensity, line-of-sight velocity, and 
longitudinal magnetic flux were obtained in the Nil 6768 



A line at a rate of one per minute for several hours. The 
linc-of-sight velocity images were available in a binned 2x2 
format, i.e. with an effective spatial scale of 1.2"/pixel. 

We observed the small Active Region (AR) NOAA 
7984 over 5 consecutive days (Aug. 15 - 19, 1996). Its ac- 
tivity was very low and although some stronger magnetic 
structures (a small spot and some pores) were present in 
the field of view (FOV), no major eruptions of magnetic 
flux were recorded. The situation was then appropriate 
to study and characterize the properties of NBPs visible 
within and around the AR. We remark that our set of 
data allows us to directly compare the NBPs as visible 
at different wavelengths with the corresponding magnetic 
field structure, and to follow them from the deep pho- 
tosphere to the higher chromosphere. In this paper we 
analyze the data obtained on August 15th, 1996, since 
for that day we had the best uniformity in time coverage 
for NSO and MDI data. Fig. 1 shows the FOV at several 
wavelengths. The period of best seeing for ground-based 
observations ran from 15:15 to 16:05 UT. This interval 
is adequate for the study of network points, since it al- 
lows an analysis of their (possible) peri odical properties , 



while they still maintain th eir identity ( |Lites et al. 1993 
von Uexkiill k Kneer 1995|) . MDI data were available for 



many hours around this interval; we consider in this work 
the period 14:00 - 17:00 UT. 

The data were re-scaled to the 0.605"/pixel of the MDI 
maps, and the co-alignment of the whole dataset was ob- 
tained comparing the position of the prominent solar fea- 
tures, i.e. the little spot and pores (see Fig. |l|). At each 
given time the alignment among the images acquired with 
different instruments was better than about 1", the mean 
spatial resolution limit of the ground-based frames. 



3. Characteristics of NBPs 

In order to identify suitable network points, we first se- 
lected on the spectra bright features showing strong K2 
peaks and enhanced wings emission in Call K. Since the 
NBPs lifetime is typically longer than 10 minutes (Rutten 
& Uitenbroek 1991), we further required that the points 
be visible for the entire observing period in the NaD2 im- 
ages as bright structures with intensity above the average. 
A total of 11 NBPs with the required characteristics were 
selected. They are distributed over the FOV both near 
the center of the AR and away from it, as can be seen 
in Fig. |l|-a. Several more network points (or structures) 
are visible on the longitudinal magnetic flux maps, but in 
this work we limit our analysis to only these 11 points for 
which corresponding Call K spectra are available. 



3.1. Intensity and velocity 

Within the areas enclosing the NBPs defined above, we 
tried to identify the bright points at each wavelength by 
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Table 1. Summary of the observations. 



Instrument 



UBF 



Zeiss 

White Light 

HSG 

MDI 



FOV 



2' X 2' 



2' X 2' 
2' X 2' 
0.75" X 2' 
10' x 6' 



Spat, resol. 
0.5" x 0.5" 



0.5" x 0.5" 
0.5" x 0.5" 
0.75" x 0.36" 
0.6" x 0.6" 



Observing A (A) 

5889.9 (NaD 2 ) 
6562.8 (Ha) 
6561.3 (Ha - 1.5 A) 
6564.3 (Ha + 1.5 A) 
5500 

3904-3941 (Call K) 
6768 (Ni I) 



FWHM (A) At (s) 



0.2 

0.25 

0.25 

0.25 

100 

0.035 

0.1 



12 



(id 



choosing an intensity (or other signature) threshold that 
could clearly separate them from their surroundings. 

On white light and Nil pseudo-continuum images, an 
intensity threshold cannot clearly discriminate between 
NBPs and other areas. The contrast averaged over the 
spatial locations corresponding to the NBPs is of the or- 
der of 1%, i.e. smaller than the rms noise calculated in 
quiet areas (about 3.5% and 2% for the white light and 
pseudo-continuum images, respectively). This is consistent 
with the observations of Topka et al. (1997), that find low 
continuum intensity contrast in network points with mag- 
netic flux density smaller than m 300 G (as typical of our 
points, see Sect. |3~2| ). 

The line of sight velocities averaged over the NBP areas 
are small, about 80 m/s downward with respect to the 
average values over the quiet areas of the FOV. Since the 
standard deviation of the measurements is about 200 m/s, 
for both NPBs and quiet areas, it is not possible to define 
a velocity threshold that allows to discriminate the NBPs 
within the FOV. The small average red-shift is however in 
agreement with recent observations by Solanki (1993) and 
Martmez-Pillet et al. (1994). 

The NBPs are instead well visible in the images ac- 
quired in the Ha far wings and NaD2 center, as sharp and 
isolated bright structures of comparable size (3"— 4" wide). 
The bright points, as seen at these wavelengths, spatially 
coincide within the overlapping error of 1". In the Ha cen- 
ter images the morphology is less clear. The presence of 
contiguous features, at times brighter than the selected 
NBPs, made more uncertain their identification. In total, 
however, 10 of the 11 considered NBPs were unambiguosly 
identified in the Ha center images. 

The characteristics of the NBPs for different signa- 
tures, including their typical contrast with respect to quiet 
areas, are summarized in Table Q It must be remarked 
that the formation height of these signatures has been 
computed in a mean quiet atmosphere, i.e. that it rep- 
resents only a generic indication for magnetic structures 
such as the NBPs. In particular, due to the Wilson depres- 
sion, the radiation coming from photospheric magnetic 
structures is believed to be formed in deeper geometric 
layers with respect to non-magnetic ones. Since the spatial 
resolution for our observations is not sufficient to resolve 
the (supposedly) elementary magnetic fluxtubes (with di- 
mensions smaller than 0.3", as seen for example in G-band 



Table 2. Characteristics of the NBPs at different at- 
mospheric heights. The first column gives the formation 
height for each wavelength, computed in the quiet Sun 



VALC model (Vernazza et al. 1981), taking into account 
the different widths of the filters. The second column 
gives the typical diameter of the structures when they 
are clearly visible. The third column gives the intensity 
contrast values, averaged over time. A range of values is 
reported when variations are large over the sample. 



Observed features 


height (km) 


size( ) 


AI/I 


White Light 







.01 


Ha wings 


100. 


3 


.04 


Pseudo-continuum (Ni I ) 


100-200. 




.01 


NaD 2 


600. 


3 - 4 


.1 


Ha center 


1500. 


4 


.1 - .2 



images), the signals we analyze are a non- linear combina- 
tion of magnetic and non-magnetic ones. 

3.2. Magnetic structures 

The selected NBPs are clearly recognizable on the 
MDI magnetic maps as sharp and isolated structures, 
3"— 4" wide. Each one corresponds to a patch of defi- 
nite polarity, with magnetic flux densities ranging from a 
minimum of 30 G, to a maximum of about 250 G (for com- 
parison, the spot in the FOV has a maximum flux density 
of 1100 G). The error on a single pixel in each image is 
given at about 15 G ( Schrijver et al. 1997). Since the nct - 
work fields are mostly vertical (see, e.g., Lites et al. 1999 ), 
the flux measure is only slightly affected by the position 
of the FOV on the solar disk (cos6» ~ 0.9). 

The magnetic evolution of the NBPs is quite varied. 
Some of the points maintain stable positions and flux val- 
ues, while others experience a steady increase during the 
observing period. In one single case we see a weak mag- 
netic structure appearing during the observing sequence, 
simultaneously with the appearance of a network bright 
point. 

The spatial correspondence between the NaD2 network 
points and the magnetic structures is very good, and will 
be analyzed in more detail in the next section. The same 
correspondence is noticeable also for the network points 
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Fig. 1. AR NOAA 7984 observed on Aug. 15, 1996, N29E16. a) NaD 2 image averaged over the time interval 15:15 - 
16:05 UT. The black squares indicate the position of the NBPs analyzed in this work, b), c) and d) Same FOV as a), 
observed in Ha, Nil pseudo-continuum and longitudinal magnetic flux. 



as seen in the Ha wings, although this is less evident than means of appropriate smoothing techniques, short term 



for NaD 2 due to their lower contrast. 



3.3. Magnetic structure and NaD 2 emission 

We checked the spatial correspondence between the net- 
work points as visible in NaD 2 and in magnetic maps. 
To this end we first removed from both signatures, by 



variations due to noise and oscillations (especially in the 
5-minutes range). We then subtracted a threshold, chosen 
as the average quiet area value for the NaD 2 images, and 
as the nominal data noise (15 G) for the magnetic ones. Fi- 
nally, within the 11 selected areas, a NBP was identified in 
both signatures as the locus of the pixels exceeding 50% 
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Fig. 2. Longitudinal magnetic flux maps (scaled between +100 and —200 G) for two different NBPs at two different 
times. The contours represent the NaD2 threshold defining a network bright point. One can see the change in shape 
and position of the magnetic structures during this interval, and the corresponding changes in NaD2 emission. The 
arrow indicates a new NBP appeared simultaneously with a weak magnetic structure. 



of the local maximum. This allowed its clear separation 
from the surroundings. 

We find that, at any given time, the NaD2 network 
points are coincident in position, size and shape with the 
corresponding magnetic patches, within 1" (the overlap- 
ping error). Any change in the characteristics of the NaD2 
NBPs reflect almost perfectly those of the magnetic fea- 
tures, within the temporal resolution of the MDI data. 
This is well exemplified in Fig. 0, where we show NaD2 
contours overlaying magnetic flux maps of several NBPs 
at two different times. To our knowledge, this is the first 
time that such a correspondence is reported at high spa- 
tial and temporal resolution. A good agreement between 
the chromospheric network emission pattern and the loca- 
tions of enhanced magnetic flux had been noted in earlier 
works (Skumanich et al 1975; schrijvcr ct al. 1989; Nin- 



dos & Zirin 1998 ) but mostly for the Call K emission, 
with lower spatial and temporal resolution. No tempo- 
ral resolution was available in the observations of Beckers 
(1976) or in those of Daras-Papamargaritis & Koutchmy 
(1983), that established a correlation between magnetic 
structures and facular structures in the wings of the Mgb 
lines. 

In analogy with the Call H and K case, we could es- 
tablish a quantitative relationship between the NaD2 ex- 
cess and the magnetic flux density for the network points. 
This property implies that also the emission in the center 
of the NaD2 line can be used as a proxy for the magnetic 
field structures. For sake of simplicity, the details of the 
determination of this relationship are given in [a|. 



4. Temporal development: light curves 

To study the temporal development of the NBPs, we com- 
puted the light curves for the bright points at each wave- 
length or signature. The curves for each NBP were ob- 
tained by selecting an area that contained the bright point 
throughout the whole observing period (even if it moved 
spatially), and then averaging, for each time, over all the 
pixels whose i nten sity exceeded the threshold value de- 
scribed in Sect. 3T. A threshold equal to the nominal data 
noi se was used for the magnetic curves. As said in Sect. 
|3.l| , the NBPs are not directly visible in some photospheric 
signatures, such as white light, Nil pseudo-continuum and 
velocity images, hence we couldn't use an intensity (veloc- 
ity) threshold to obtain the corresponding light curves. To 
guarantee the comparability with the other light curves, in 
these cases we computed, at each time, the average value 
over the spatially corresponding magnetic areas. 

We also computed the light curves of 11 areas ran- 
domly selected in the quiet regions of our FOV, and of 
size comparable to that of the NBPs (about 3"x3"). No 
threshold was applied for their computation. These quiet 
areas should represent the so-called internetwork regions, 
which appear field-free at MDI sensitivity. Using the same 
number of internetwork areas as of network bright points, 
and performing the analysis in the same way, will give 
us confidence in the comparison in a statistical sense. We 
remark that this is not always the case, especially for spec- 
trographical observations where the slit samples a number 
of internetwork points that is usually much larger than 
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Fig. 3. Temporal evolution for a single NBP of velocity, 
magnetic flux density, intensity in Ha wings, NaD2, Ha. 
Long-period trends have been subtracted using a smooth- 
ing window of 600 s (see next section) . Intensity variations 
are after subtracting a value of 10 4 arbitrary units 

that of network structures. In Fig. ^, as an example, we 
show the light curves of different signatures for a NBP. 
Fluctuations are evident at each wavelength. The red and 
blue wings of Ha show a very similar temporal evolution, 
with simultaneous intensity variations of the same mag- 
nitude for most NBPs. For other signatures a direct com- 
parison does not show any simple relation between the 
variations at different atmospheric heights. 

We must comment here on the fact that most authors, 
when studying the temporal properties of network points, 
use a different method to determine the "light curves". 
Basically, the network points are spatially identified using 
the temporal average of some suitable signature, and then 
the temporal development of each image pixel belonging 
to the average structures is considered and analyzed. We 
believe that the method we adopted has several advan- 
tages: - Since the network points might move over the 
course of time, they do not always correspond spatially to 
the structures identified on the average maps. Our method 
guarantees that the structure is properly followed in time, 
avoiding the loss of relevant pixels, or the inclusion of spu- 
rious ones; - If the magnetic structure giving rise to the 
network point inclines with height, a set of pixels identify- 
ing the NBP at a given wavelength might not well repre- 
sent the same structure at a different wavelength. This is 
especially relevant when performing comparisons between 
different atmospheric layers, for example in the phase dif- 
ference analysis of Sect. 0. We overcame this problem by 
selecting areas large enough to include the NBPs at each 
wavelength and each time, as explained earlier. 

Computing the light curves as an average over a given 
area implies the assumption of a spatial coherence over 
the whole area (of about 3"x3", equivalent to a spatial 



frequency of about 3 Mm -1 ) competing to each NBP. This 
seems a reasonable assumption because we do not see any 
significant inhomogeneities within the single structures. 

5. Power spectra 

A search for possible periodicities in the fluctuations of the 
NBPs light curves was performed using temporal power 
spectra. Before computing the power spectra, the light 
curves were detrended using a smoothing window of 600 s. 
A check on this procedure showed that changing the smo- 
othing window between 360 s and 840 s affected the power 
at frequencies lower than 1.2 mHz, but without changing 
the frequency of the peaks. The power at higher frequen- 
cies remained unaffected. 

A power spectrum was computed for each light curve 
of the 11 NBPs and of the 11 internetwork areas. To an- 
alyze the differences between these two atmospheric com- 
ponents, we averaged separately the power spectra over all 
of the NBPs and over all of the quiet regions. In Fig. ^ we 
show some of these averaged curves. It must be remem- 
bered that the NSO and MDI observations have different 
temporal coverage (50 and 180 minutes, respectively) and 
temporal resolution (12 s and 60 s), so that lower tem- 
poral frequencies are better represented in the Nil series, 
while frequency coverage extends to higher frequencies for 
ground-based data. However, the power at v > 8 mHz in 
all the ground-based signatures is due essentially to noise 
(Fig. |J), hence the analysis can be limited to the frequency 
coverage of the Nil observations, 0-8 mHz. 

We describe here the power spectra characteristics 
from lower photospheric signatures to higher chromo- 
spheric ones. Intensity fluctuations may be plausibly inter- 
preted as temperature fluctuations for photospheric LTE 
signatures such as Ni I or the Ha wings, formed over 
depths where the velocity gradient is small. In these cases, 
the intensity fluctuations directly reflect fluctuations of 
the source function (the Planck function) and hence of 
temperature. In the center of chromospheric lines as NaZ?2 
or Ha the NLTE effects are important and the intensity 
fluctuations are the response to temperature, density and 
even velocity changes (Cram 1978) and then are a sort 
of average of the variation of the state variables of the 
chromosphere. 

5.1. Photospheric signatures 

In Fig. | a, b, d we show the power spectra of Ni I pseudo- 
continuum intensity, Ni I velocity and Ha red wing inten- 
sity. The power spectra computed for white light intensity 
(not shown in figure) and for the Nil pseudo-continuum 
are very similar, and we are confident that we can com- 
pare the two series of observations, even if their frequency 
resolution is different. The Ha blue wing spectrum (not 
shown in figure) has a similar behaviour but a lower am- 
plitude with respect to the red wing (see Table ||) . A value 
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Fig. 4. Power spectra in arbitrary units, averaged for NBPs (solid) and internetwork regions (dashed). The 50% 
confidence limit is given by the range 0.8— 1.5x the actual power values; for clarity error bars are not plotted in the 
figure. 



smaller than 1 for the ratio of the power in the blue and 
red wings of Ha is consistent with the observations of 
Bcrtcllo (1987), that found the same trend for the power 
of velocity oscillations in the wings of photospheric lines 
formed at heights lower than 150 km. 

As is well known, the distribution of power for the pho- 
tospheric velocity fluctuations is rather different from the 
one of pseudo-continuum intensity fluctuations (Fig. ^ a- 
b). The velocity power is concentrated around the range 
of frequencies corresponding to the 5-minutes oscillations. 
The pseudo-continuum power spectrum peaks at low fre- 
quencies, around 1.5 mHz and then show a decay that 
might indicate the stochastic character of the granulation 
intensity variation, as already reported for the first time 
in Noyes (1967). 



As a global characteristic, power spectra computed in 
photospheric signatures do not show any significant differ- 
ence between network and internetwork structures within 
the 50% confidence limit (Fig. [| a, b, d). This result is 
consistent with previous spectral observations by several 
authors ( Dcubner fc Fleck 1990] ; Kulaczewski 1992 ; Lites 
et al. 1993) that analyzed both intensity and velocity os- 
cillations in the photosphere for network and internetwork 
features. 

The power spectrum of the magnetic flux variations 
averaged over the NBPs is shown in Fig. ^ c. Internet- 
work areas are not considered because the noise in the 
magnetic flux measure is too high for a reliable determi- 
nation of fluctuations. Significant peaks are visible at low 
frequency (around 1.5 mHz), indicating long term evo- 
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lution of the magnetic field, and around 3.5 mHz corre- 
sponding to the 5 minutes oscillations. A signal at the 
latter timescale might represent the magnetic response to 
oscillations already present in the photosphere and be of 
importance in the context of generation and dissipation of 



Table 3. Power values in arbitrary units for various ob- 
served features. A range of values is reported when varia- 
tions are large. 



MHD waves in the solar atmosphere (Ulrich 1996). Obser- 
vations of flux variations in small magnetic structures are 
scarce in the literature, but we can compare this result 
with those presented by Norton et al. (1999), that used 
a similar set of MDI data obtained in the area of a big 
sunspot. They found a significant peak near 5 min only 
for structures whose magnetic flux density exceeded 600 
G, while the points we analyzed had a maximum value of 
about 300 G. 



5.2. Chromospheric signatures 

The intensity power spectra computed in chromospheric 
signatures display strong differences between NBPs and 
internetwork as shown in Fig. ^ e - f . We will analyze in 
detail these differences, keeping in mind that the regime 
of oscillations changes with height in the chromosphere. 

First of all we notice that the power spectrum of inter- 
network intensity fluctuations in NaD2 shares some char- 
acteristics with the photosphcric Nil velocity power spec- 
trum rather than with the one of Ni I intensity. In par- 
ticular the strongest peak appears around 3.5 mHz, while 
the enhanced low frequency component, typical of photo- 
spheric intensity power spectra, is lacking. This charac- 
teristic could be explained if the intensity fluctuations in 
the NaD2 line center were related more to velocity than 
to temperature perturbations. This might be indirectly 
confirmed by the results of Pallc et al. (1999) in their 
study of the current performances of the GOLF experi- 
ment on SOHO. In determining the relative contributions 
of velocity and intensity signals to the intensity variations 
measured in the blue wings (—100 mA) of the sodium dou- 
blet, they conclude that the effect due to "pure" intensity 
changes is only 14% that of velocity changes for the p - 
mode frequency range. Since the width of the filter used 
for our observations includes that same portion of the line 
wing, this conclusion might apply, at least partially, also 
to our case. 

Comparing network and internetwork power spectra 
for NaD2 , we see that the power of the NBPs is smaller 
than the corresponding power for internetwork at each 
temporal frequency. In particular, even if both NBPs and 
internetwork points show a maximum around 3.5 mHz, 
the power at this frequency is almost an order of mag- 
nitude smaller in the NBPs. (Fig. || d suggests that this 
effect might be already present in the wings of Ha, al- 
though its amplitude is not large enough to give an un- 
ambiguous result). This suppression of power in the NBPs 
indicates that the presence of the magnetic field in some 
way perturbs and reduces the oscillations at low chromo- 
spheric levels, especially in the p-mode range. A compres- 
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sion of power in magnetic structures at frequencies below 
7-8 mHz, for lines formed at similar heights, has not been 
reported by other authors. Only Al et al. (1998) observe a 
similar effect, but much smaller, for the power spectrum of 
velocity fluctuations measured in the center of NaD2 with 
a narrowband filter (30 mA). We think that the stronger 
effect seen in our observations is real because our anal- 
ysis selects the horizontal scale typical of chromospheric 
network to determine the light curves and the power spec- 
tra (see Sect. [|), and is therefore more suitable to outline 
characteristics and differences on the same horizontal scale 
for NBPs and internetwork areas. 

The power spectrum computed for the intensity of Ha 
center is shown in Fig. [| f. For both NBPs and inter- 
network the power distribution peaks at low frequencies 
{v < 3 mHz), without any relevant peak at the p-mode 
frequencies. No enhanced power for "3-minutes" oscilla- 
tions [y > 5 mHz), is detectable in the internetwork. This 
is consistent with observations in the Ha center by Cram 
(1978) and in the Ca II - H3 by |Lites et al. 1993| showing 



a power peak in the "3 minutes" range only for velocity 
power spectrum. 

The spectrum of NBPs in Ha line has a power higher 
than that of the internetwork areas at each frequency, re- 
versing the effect present in the NaD2 line, and shows 
three well separated peaks reminiscent of the peaks ob- 
served by Lites et al. (1993). The more relevant peak is 
around 2.2 mHz ("7- minutes" oscillations) and is lowered 
about a factor 6 in the internetwork. We cannot judge 
on the relevance of the peak around 1.3 mHz, since its 
amplitude is heavily affected by the smoothing window, 
as described in SectJ^, and we cannot consider real the 
peak at 0.6 mHz, because it is related to the time interval 
of our observations. However, the increasing power at low 
frequencies in the spectrum of NBPs suggests that the role 
of the magnetic field in the oscillations, detectable at high 
chromospheric levels, is certainly different from the one at 
lower levels and strenghtens the hypothesis of magnetic- 
hydrodynamic waves present at these high levels. 

Characteristics of power spectra for network and inter- 
network are summarized in Table [| for the two relevant 
frequency windows 1.5 - 2.5 mHz and 3-4 mHz. 
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6. Phase difference and coherence spectra 

In order to look for propagation characteristics of waves 
at different heights in the atmosphere, we computed phase 
difference (A<£>) and phase coherence (C) spectra for many 
signature pairs (Straus 1995). We exclude in this analysis 
the NaD2 signature because, if the measured intensity os- 
cillations are due essentially to velocity oscillations, we 
cannot establish the direction of the motion and hence 
assign the correct value to the phase difference. 

Following Edmonds & Webb (1972), we adopted a 
smoothing width of about 1.5 mHz in the Fourier do- 
main for the computation of phase and coherence spectra. 
Frequencies smaller than 0.7 mHz hence do not convey 
any significant information. Taking into account our smo- 
othing width, a coherence smaller than ^0.5 implies that 
phase differences at those frequencies are completely un- 
reliable. 

NSO and MDI data were analyzed independently, us- 
ing their own frequency resolution and coverage. Finally, 
the study was performed separately for internetwork and 
network areas, to distinguish between features with differ- 
ent magnetic characteristics. We describe the characteris- 
tics of the spectra from lower photospheric layers through 
higher chromospheric ones. 

In Fig. ^ left and central column, we show intensity 
(I— I) phase difference and coherence spectra for the pairs 
Ha + 1.5 A /Ha - 1.5 A, and Ha + 1.5 A / white light, 
originating at photospheric levels. We use here the white 
light signal (and not the Nil pseudo-continuum) in order 
to keep the maximum possible frequency resolution. Phase 
difference and coherence spectra between magnetic flux 
density and velocity (B— V) for the NBPs are shown in 
Fig. right column. 

To search for the relationship between the oscillations 
present at chromospheric and photospheric levels we com- 
puted the I— I phase difference and coherence spectra sep- 
arately for the pairs Ha center / Ha —1.5 A and Ha center 
/ Ha +1.5 A, shown in Fig. |[ At each considered atmo- 
spheric level, a general characteristic is that the coherence 
for NBPs is smaller than for internetwork, hence the phase 
values for NBPs are more uncertain. 

We examine different signature pairs separately in the 
two frequeny intervals 1.5 - 2.5 mHz and 3-4 mHz, dis- 
regarding the features with negligible power, and we sum- 
marize in Table || the phase and coherence values for each 
pair. 



6.1. Low frequency (1.5 - 2.5 mHz) 

For both internetwork and network points the I— I phase 
difference between the Ha red and blue wings is 0°. The 
two signals should be in phase if the observed oscillations 
are due only to temperature and in antiphase if due to 
velocity. We can then state that the observed oscillations 
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at low frequencies are essentially due to temperature os- 
cillations. 

It follows from the previous considerations that the 
analysis of the phase difference spectra between the 
Ha red wing and white light is essentially a study of the 
correlation between temperature oscillations at slightly 
different levels (Ah < 100 km in the quiet average photo- 
sphere). For internetwork areas the extremely high coher- 
ence makes the A<!> value (5 - 10°) highly significant (see 
Table ||) and strongly suggests that the observed power is 
due to oscillations. A positive value of A<!> between two 
layers with decreasing heights, indicates the presence of 
waves directed radially inward. For the internetwork ar- 
eas, free from magnetic fields, at the spatial frequency of 
about 3Mm _1 used in our analysis and within the con- 
sidered frequency range, these waves might be interpreted 
as gravity waves (Straus 1995). This same interpretation 
has been adopted for internetwork by Rutten (2000), who 
reported a similar value of A<& between the intensity of 
two continuum levels observed by TRACE. Downward di- 
rected gravity waves had been proposed earlier by Staiger 
et al. (1984) to explain similar phase differences for veloc- 
ity signatures at photospheric levels. 

At chromospheric levels, the most significant peak in 
the power spectrum of the Ha center intensity appears 
at 2.2 mHz for NBPs and is not related to the oscilla- 
tions present at photospheric levels ( C = 0.45 between 
Ha center /Ha +1.5 A, see Table ||). This means that at 
chromospheric levels the NBPs experience a new regime of 
oscillations that seem to be independent from what hap- 
pens in photosphere. 

6.2. p-mode frequency (3.0 - 4-0 mHz) 

In the p-mode range of frequencies, for internetwork there 
is a 10°phase lag between the two Ha wings (see Table f§ 
that can be due to different coupling of velocity and in- 
tensity fluctuations in the two wings. This effect has been 
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Fig. 5. Top and mid panels show the phase difference I — I, as a function of frequency, of photosphcric signatures 
averaged for network and internetwork regions. Vertical bars, representing la error associated with the mean phase 
value, are clipped if higher than the y-axis limits. The bottom panels show the phase coherence spectra (computed 
with a smoothing width of 1.5 mHz, see text) averaged for network (solid line) and internetwork (dashed line) regions. 



extensively studied for photospheric lines by Cavallini et 
al. (1987) and Alamanni et al. (1990). 



As described in Sect. |5^lJ, the presence of a peak at 
these frequencies in both the magnetic flux and velocity 
power spectra could betray the presence of MHD waves 
in the network structures. In the limit of ideal MHD, a 
definite phase relation between velocity and magnetic field 
variations is expected as signature of Alfven waves ( A<I> = 
0°) or magnetoacoustic waves (A$ = 90°, with v leading 
B. See Ulrich 1996). However, we find a very low coherence 
(< 0.4) at all frequencies in the NBPs, i.e. the magnetic 
fluctuations are not related to the velocity ones, at least 
with the present sensitivity and resolution. 



7. Discussion and conclusions 

The observations presented in this paper allowed us to de- 
fine the characteristics of network bright points at differ- 
ent atmospheric heights, and to compare them with those 
of the surrounding internetwork areas. We improved on 
the existing statistics using a good-sized sample of NBPs, 



and the same number of "test" internetwork areas, de- 
fined in a comparable way. The method we adopted to 
study the temporal evolution of NBPs insures that each 
bright structure is properly followed in time and position 
at each height. In fact, the evaluation of the light curves 
and their properties after a spatial averaging over a well 
defined area guarantees that we are studying the same 
NBP at all heights, and avoids the problem (first pointed 
out by Lites 1994) of a possible structure displacement 
due to the magnetic field inclination. Given the charac- 
teristic horizontal size of the NBPs, the analysis and the 
comparison of power spectra and phase differences con- 
cern the propagation of waves pertaining to a horizontal 
wavenumber of about 3 Mm . 

The quasi-simultaneous series of NaD2 images and of 
MDI maps allowed us to establish for the first time a cor- 
respondence between NaD2 bright network and magnetic 
network at high spatial and temporal resolution. A corre- 
spondence between bright chromospheric structures (Ca 
II, Lya, Mg I and UV continuum) and magnetic struc- 
tures had been observed before, but not at this high tern- 
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Fig. 6. Top and mid panels show the phase difference 
I — I, as a function of frequency, for the pairs Ha cen- 
ter / Ha —1.5 A and Ha center / Ha +1.5 A averaged 
for network and internetwork regions. Vertical bars, rep- 
resenting la error associated with the mean phase value, 
are clipped when higher than the y-axis limits. The bot- 
tom panels show phase coherence spectra (computed with 
a smoothing width of 1.5 mHz, see text) averaged for net- 
work (solid line) and internetwork (dashed line) regions. 



poral resolution. We also established for the NBPs a quan- 
titative relationship between the Na excess and the cor- 
responding absolute value of magnetic flux density. This 
relationship is best expressed by a power law with an expo- 
nent very close to the one found by Schrijver et al. (1989, 
1996) for the Ca II - K excess, and indicates that the emis- 
sion in NaD2 may be used as a proxy for the magnetic flux 
density. 

The NBPs considered in this work have the follow- 
ing properties: - are bright in the Ca II wings and in 
the Ca II K 2 peaks; - are visible in the NaD 2 images 
for about 1 hr; - coincide spatially with the magnetic 
structures; - are nearby or within a lower activity region. 
The general characteristics found for these NBPs do not 
differ from the ones derived in absolutely quiet regions 
( [Deubner fc Fleck 1990| , |Lites et al. 1993[ ). 

Our results referring to photospheric and chromo- 
spheric properties are so summarized: 

At 'photo spheric levels: No difference is detected be- 
tween network and internetwork power spectra, either in 



intensity or in velocity, within the limits of sensitivity and 
accuracy of the instruments used for this work. The phase 
difference spectra between photospheric signatures in gen- 
eral do not show different characteristics for network or in- 
ternetwork. However, when analyzing the phase difference 
between Ha red wing and white light images(Aft, < 100 
km), we find 5° < A$ < 10° in the frequency window 
1.5 - 2.5 mHz and in the internetwork. (The A(f> value 
is more uncertain in the network, due to a lower coher- 
ence value). A phase lag of this amplitude and sign is 
usually considered a signature of gravity waves directed 
radially inward. A possible explanation for their origin 
might be sought in recent models of convection, described 
as a non-local process driven by cooling at the solar sur- 
face rather than by heating from the lower layers (Spruit, 
1997). One can imagine that the downward flowing cooled 
plasma can trigger some inward directed waves, and hence 
justify the fact that the external layers "lead" the deeper 
layers. The general inhibition of convection in magnetic 
structures might be the reason for the lack of this signa- 
ture in network points. 

The power spectrum of the magnetic flux variations in 
NBPs shows a small but significant peak around 3 mHz, 
that could be related to a "transformation" of acoustic 
waves into MHD waves. However, the phase difference 
and coherence spectra between magnetic flux and velocity 
(B— V) for the NBPs indicate a very low correlation be- 
tween the two signals so we cannot conclude anything on 
the presence of MHD waves within the network points. 

At chromo spheric levels: Network and internetwork ar- 
eas have a rather different behaviour in the power spectra. 
We do not see any evidence for the typical chromospheric 
period of 3 minutes (but it must be reminded that they 
are best seen in velocity variations rather than intensity) . 
In the low chromospheric levels, where NaD 2 originates, 
the NBPs power spectrum is compressed at all frequencies 
if compared to the internetwork, while in the high chromo- 
sphere, where Ha originates, the power of NBPs is higher 
than the one of internetwork. This opposite effect may be 
an indication that the magnetic field disturbs and reduces 
the amplitude of oscillations already present in the low 
chromosphere while it assumes a leading role in the high 
chromosphere. 

In the layers contributing to the NaD 2 emission it 
seems that the oscillations present in network points 
change regime with respect to both the photosphere and 
the high chromosphere and we think that it would be im- 
portant to perform observations of NBPs in tha Na line, 
with high spectral resolution. Unfortunately we cannot 
analyze the phase difference spectrum for NaD 2 intensity 
fluctuations with respect to others formed at different lay- 
ers, since the NaD 2 intensity fluctuations, measured with 
the UBF filter (FWHM= 0.2A), are more related to ve- 
locity than to temperature fluctuations (see Sect |5.2j ) . 

The power spectrum of Ha intensity in NBPs has the 
more relevant peak at 2.2 mHz, but this signal is not cor- 
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related with the photospheric fluctuations, as indicated 
by the very low coherence measured at all frequencies be- 
tween the Ha core and the blue and red wings. We can 
then confirm, using a larger sample of NBPs, the presence 
of the peak found by Lites et al. (1993) around 2 mHz 
in the power spectrum of K3 velocity fluctuations for one 
network point. Kalkofen (1997) and Hasan & Kalkofcn 
(1999) proposed an explanation for this peak in terms of 
transverse magneto-acoustic waves in magnetic flux tubes, 
excited by granular buffeting in the solar photosphere. In 
their model the low coherence between photospheric and 
chromospheric signatures could be explained by a partial 
conversion of the transverse waves to longitudinal modes 
in the higher chromosphere. 

A general result of our analysis, valid from the low pho- 
tosphere to the high chromosphere, is that the NBPs al- 
ways show a coherence lower than the internetwork, point- 
ing out that the presence of the magnetic field changes 
the propagation regime of waves with respect to the non- 
magnetic regions. 
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Appendix A: A relationship between Na excess 
and magnetic field 

A quantitative relationship between the Call K line 
core intensity and the absolute value of the mag- 
netic flux density has been clearly established by sev- 



eral studies (Skumanich et al 1975; Schrijver et al. 1989; 
Nindos & Zirin 1998). The best agreement with the data 



is given by a power law relation with an exponent of 0.6 
( [Schrijver et al. 198E ). Given the excellent spatial coin- 
cidence of the sodium network points and the magnetic 



structures present in our data (see Sect. 3.3, we searched 
for such a quantitative relationship also for this case. 
Fig. A.l shows a scatterplot of the sodium "excess" 



for the 11 NBPs, vs. the corresponding absolute values of 
magnetic flux density. As done by Schrijver et al. (1989, 
1996) for the Call K emission, we subtracted to the NaD2 
intensities a threshold equal to the average intensity in the 
quiet areas (10 in arbitrary units). The graph was ob- 
tained plotting the temporal average of both the sodium 
intensity and magnetic flux. It does not display, hence, 
any temporal variation of either quantity, but only a gen- 
eral trend among the persistent structures. Five-minutes 
oscillations do not play a role in this relationship, nor do 
they contribute to the scatter of the points. Saturation, as 



1000 - 



X 



0-! 
Q 




100 i 



10 100 
Magnetic Flux Density (Gauss) 

Fig. A.l. Scatter diagram of the NaD2 intensity, after 
subtraction of a threshold of about 10 4 (see text), vs. the 
absolute longitudinal magnetic flux density. All the 11 an- 
alyzed NBPs are plotted. Solid line is the power law fitting 
described in the text. 

reported by Schrijver et al. (1989) is not apparent in our 
data, but it must be remarked that the maximum value 
of the magnetic flux density for the average quantities is 
below 400 G, the "critical" value indicated by those au- 
thors. 

The sodium excess is best fitted by a power law of 
the type Al cx $0, with /3=0.58 ±0.1. The exponent is 
very close to /3=0.6±0.1 found by Schrijver et al. (1989, 
1996) for the Call K excess. These results indicate that 
the emission in the center of the NaD2 line is also a good 
proxy for the magnetic flux density and, at least for values 
of magnetic flux density smaller than a few hundred G, its 
use is equivalent to that of the Call K emission. 
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